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ABSTRACT 

This is the first of a series of papers devoted to the investigation of a large sample of 
brightest cluster galaxies (BCGs), their kinematic and stellar population properties, 
and the relationships between those and the properties of the cluster. We have obtained 
high signal-to-noise ratio, long-slit spectra of these galaxies with Gemini and WHT 
with the primary purpose of investigating their stellar population properties. This 
paper describes the selection methods and criteria used to compile a new sample 
of galaxies, concentrating on BCGs previously classified as containing a halo (cD 
galaxies), together with the observations and data reduction. Here, we present the full 
sample of galaxies, and the measurement and interpretation of the radial velocity and 
velocity dispersion profiles of 41 BCGs. We find clear rotation curves for a number 
of these giant galaxies. In particular, we find rapid rotation (> 100 km s -1 ) for two 
BCGs, NGC6034 and NGC7768, indicating that it is unlikely that they formed through 
dissipationless mergers. Velocity substructure in the form of kinematically decoupled 
cores is detected in 12 galaxies, and we find five galaxies with velocity dispersion 
increasing with radius. The amount of rotation, the velocity substructure and the 
position of BCGs on the anisotropy-luminosity diagram are very similar to those of 
"ordinary" giant ellipticals in high density environments. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: kinematics and dynamics 
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1 INTRODUCTION 

The galaxies in the centres of clusters are unique. They are 
usually the dominant, brightest and most massive, galaxies 
in their clusters. A wealth of imaging data has been accu- 
mulated for these intriguing objects (Malumuth & Kirshner 
1985; Schombert 1986, 1987, 1988; Postman & Lauer 1995; 
Collins & Mann 1998; Brough et al. 2002; Laine et al. 2003), 
however the spectroscopic data is limited to very small sam- 
ples or narrow wavelength coverage. 

Tonry (1984), Tonry (1985), Gorgas, Efstathiou & 
Aragon-Salamanca (1990), Fisher, Illingworth & Franx 
(1995a), Fisher, Franx & Illingworth (1995b), Cardiel, Gor- 
gas & Aragon-Salamanca (1998) and Carter, Bridges & Hau 
(1999) each investigated 18 or less galaxies, and moreover 
measured only three or less indices. Brough et al. (2007) used 
a wide wavelength range but included only three brightest 
cluster galaxies (BCG). Von der Linden et al. (2007) studied 
BCGs from the Sloan Digital Sky Survey (SDSS), however 
the spatial information is lacking since fibers were used. 

Among BCGs there exists a special class of galaxies cat- 
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alogued as cD galaxies, although there exists a great deal of 
confusion over the exact meaning of the classifications gE, 
D and cD. Matthews, Morgan & Schmidt (1964) outlined 
the following definition: "D galaxies have an elliptical-like 
nucleus surrounded by an extensive envelope. The super- 
giant D galaxies observed near the centre of a number of 
Abell's rich clusters have diameters 3-4 times as great 
as the ordinary lenticulars in the same clusters. These very 
large D galaxies observed in clusters are given the prefix "c" , 
in a manner similar to the notation for supergiant stars in 
stellar spectroscopy". Later, Schombert (1987) defined gE 
galaxies as distinct from other early-type galaxies by their 
large size; D galaxies as being gE galaxies with a shallow 
surface brightness profile slope; and cD galaxies as D galax- 
ies with large extended stellar haloes. The last are also more 
diffuse than normal ellipticals (Schombert 1986) . Because of 
the confusion surrounding the definition of D galaxies, and 
the fact that they are rarely regarded as a separate type of 
object in the modern literature, this study will only refer to 
cD and non-cD BCGs (i.e BCGs containing a halo or not). 
For BCGs, we adopt the definition to comply with recent 
literature (for example Von der Linden et al. 2007), where 
BCG refers to the central, dominant galaxy in a cluster. For 
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a small fraction of clusters, the BCG might not strictly be 
the brightest galaxy in the cluster. 

Approximately 20 per cent of rich clusters of galax- 
ies contain a dominant central cD galaxy (Dressier 1984; 
Oegerle & Hill 2001), although they can be found in poor 
clusters as well (Giacintucci et al. 2007). Some clusters 
have more than one cD galaxy, but a cD galaxy is always 
the dominant member of a local subcluster. The surface- 
brightness profiles of cD galaxies are displaced above the 
de Vaucouleurs law (de Vaucouleurs 1948) at large radii. 
The break in the cD galaxy surface-brightness profile typi- 
cally occurs between 24 and 26 mag/arcsec 2 in the V-band 
(Sarazin 1988). The interpretation of this deviation is that 
the galaxy is embedded in an extensive luminous stellar halo. 
Three main theories have been proposed over the last four 
decades to explain the properties of BCG, and in particular 
cD, galaxies. 

Theory 1. In the first theory, BCG formation is caused by 
the presence of cooling flows in clusters of galaxies (Cowie & 
Binney 1977). Cooling flow clusters are common in the local 
universe and BCGs are most often found at the centres of 
these systems (Edwards et al. 2007). If the central cluster 
density is high enough, intracluster gas can condense and 
form stars at the bottom of the potential well. Observations 
that support this idea are blue- and UV-colour excesses ob- 
served in the central galaxy of Abell 1795 (indicative of star 
formation) by McNamara et al. (1996) and molecular gas 
detected in ten out of 32 central cluster galaxies by Salome 
& Combes (2003). Cardiel et al. (1998) obtained radial gra- 
dients for the D4000 and Mg2 spectral features in 11 central 
cluster galaxies. Their observations were consistent with an 
evolutionary sequence in which radio-triggered star forma- 
tion bursts take place several times during the lifetime of 
the cooling flow in the centre of the cluster. However, McNa- 
mara & O'Connell (1992) find only small colour anomalies 
with small amplitudes, implying star formation rates that 
account for at most a few percent of the material that is 
cooling and accreting onto the central galaxy. 

More recently, XMM-Newton observations showed that 
the X-ray gas in cluster centres does not cool significantly 
below a threshold temperature of kT ~ 1 — 2 keV (Jordan et 
al. 2004, and references therein). The central cluster galax- 
ies often host radio-loud AGN which may account for the 
necessary heating to counteract radiative cooling (von der 
Linden et al. 2007). Although BCGs are probably not com- 
pletely formed in cooling flows, the flows play an important 
role in regulating the rate at which gas cools at the centres 
of groups and clusters. 

Theory 2. The second theory was proposed by Merritt 
(1983) and suggests that the essential properties of BCGs, 
and in particular those with haloes, are determined when the 
clusters collapse (primordial origin). Thereafter, frequent 
mergers of galaxies would be inhibited by the relatively high 
velocities between galaxies. Merritt (1983) argued that all 
galaxies had large haloes early in the life of the cluster. 
These haloes were then removed by the mean cluster tidal 
field during the initial collapse and returned to the cluster 
potential, except for the central member which remained 
unaffected because of its special position with respect to the 
cluster potential. 



Theory 3. The third, and most widely accepted, theory is 
in the context of the ACDM cosmology and relates the for- 
mation of the central galaxy to mergings with or captures 
of less massive galaxies, and is known as "galactic cannibal- 
ism". It was first proposed by Ostriker & Tremaine (1975) 
and developed by Ostriker & Hausman (1977). 

The most complete quantitative prediction of the for- 
mation of BCGs in the now standard CDM model of struc- 
ture formation is by De Lucia & Blaizot (2007). They used 
N-body and semi-analytic techniques to study the forma- 
tion and evolution of BCGs and found that, in a model 
where cooling flows are suppressed at late times by AGN 
activity, the stars of BCGs are formed very early (50 per 
cent at 2 ~ 5 and 80 per cent at 2 ~ 3) and in many small 
galaxies. They also found that BCGs assemble late: half of 
their final mass is typically locked up in a single galaxy after 
2 ~ 0.5 (illustrated in their figure 9). A very similar conclu- 
sion was reached by Romeo et al. (2008), who performed 
N-body and hydrodynamical simulations of the formation 
and evolution of galaxy groups and clusters in a ACDM cos- 
mology to follow the build-up of two clusters and 12 groups. 
Observationally, Aragon-Salamanca et al. (1998) examined 
the K-band Hubble diagram for BCGs up to a redshift of 2 
= 1. They found that the BCGs had grown by a factor of 
two to four since z = 1. Brough et al. (2002) found a similar 
result but discovered that the mass growth depended on the 
X-ray luminosity of the host cluster. They found that BCGs 
in high X-ray luminosity clusters showed no mass accretion 
since 2 = 1 as opposed to BCGs in low X-ray luminosity 
clusters which grew by a factor of four. However, the recent 
near-infrared photometric study of 42 BCGs by Whiley et al. 
(2008) in the 0.2 < 2 < 1 range contradicts this. They stud- 
ied the colour and rest-frame K-band luminosity evolution 
of BCGs and found it to be in good agreement with popu- 
lation synthesis models of stellar populations which formed 
at 2 ~ 2 and evolved passively thereafter. 

Using the Millennium Simulation, De Lucia et al. (2006) 
studied how formation histories, ages and metallicities of el- 
liptical galaxies depend on environment and on stellar mass. 
Their Figure 9 shows the effective number of progenitors of 
early-type galaxies as a function of galaxy stellar mass. The 
number of effective progenitors is less than two for galaxies 
up to stellar masses of ~ 10 11 Mq. This function suddenly 
increases up to the value of approximately five effective pro- 
genitors for the mass of a typical BCG. 

A related theory, called tidal stripping, was first pro- 
posed by Gallagher & Ostriker (1972). Cluster galaxies that 
pass near the gravitational centre of the cluster may be 
stripped of some of their material by the tidal forces from 
the cluster potential or the central galaxy potential. The 
stripped material falls to the centre of the potential well, and 
could contribute to the observed haloes of cD galaxies. The 
most massive galaxies surrounding the central galaxy would 
be preferentially depleted as they are most strongly affected 
by dynamical friction (Jordan et al. 2004). The difference be- 
tween stripping and primordial origin is that stripping (and 
cD halo formation) begins after cluster collapse whereas pri- 
mordial origin assumes that the tidal events occur before 
collapse, and that the cD halo is not a consequence of tidal 
stripping (Schombert 1988). 

The observation of multiple nuclei in central galaxies 
favours the cannibalism theory (Postman & Lauer 1995). 
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Tonry (1984) observed the velocity and velocity dispersion 
profiles of NGC6166 and NGC7720 and their multiple nu- 
clei. He found that the stellar velocity dispersion of the cen- 
tral galaxy demonstrated that the multiple nuclei are not 
following circular orbits. This was followed by a bigger sam- 
ple of 14 multiple nuclei BCGs for which the redshifts and 
stellar velocity dispersions are presented in Tonry (1985). 
Yamada et al. (2002) showed that the BCG in a cluster 
at z = 1.26 is composed of two distinct sub-units that are 
likely to fully merge on a time-scale of 10 8 years. Jordan et 
al. (2004) studied the globular cluster systems in BCGs with 
confirmed haloes from Hubble Space Telescope observations. 
They concluded that the observed globular cluster metallic- 
ity distributions are consistent with those expected if BCGs 
galaxies form through cannibalism of numerous galaxies and 
protogalactic fragments that formed their stars and globular 
clusters before capture and disruption, although they state 
that the cannibalism scenario is not the only possible mecha- 
nism to explain these observations. The Jordan et al. (2004) 
globular cluster data also suggests that BCGs experienced 
their mergers prior to cluster virialisation, yet the presence 
of tidal streams suggest otherwise (Seigar, Graham & Jerjen 
2007). 

Carter & Metcalfe (1980) and West (1989) showed that 
the major axis of BCGs tends to be aligned with the ma- 
jor axis of the cluster galaxy distribution. Recent studies of 
the Coma cluster (Torlina, De Propris & West 2007) show 
strong evidence that there are no other large-scale galaxy 
alignments other than for the BCGs. This also confirms that 
BCGs form via mechanism related to collimated infall of 
galaxies along the filaments and the growth of the cluster 
from the surrounding large scale structure (Boylan-Kolchin, 
Ma & Quataert 2006). 

In the CDM cosmology it is now understood that local 
massive cluster galaxies assemble late through the merging 
of smaller systems. In this picture, cooling flows are the main 
fuel for galaxy mass-growth at high redshift. This source is 
removed only at low redshifts in group or cluster environ- 
ments, due to AGN feedback (De Lucia & Blaizot 2007). 

Some of the outstanding issues are: whether BCGs have 
a different formation mechanism than elliptical galaxies; and 
if the formation of BCGs are controlled by environment. To 
address these issues the present project was initiated, and 
long-slit spectra of a large and statistically significant sam- 
ple of BCGs were obtained. These data provide a set of 
spectral indices, covering a wide wavelength range. The lu- 
minous central galaxies will be contrasted with other early- 
type galaxies to look for relative differences in their evolu- 
tion, using features sensitive to stellar population age and 
abundances of various elements. 

The first part of this project entails the determination 
of the stellar kinematics through derivation of velocity and 
velocity dispersion profiles. A forthcoming paper will be de- 
voted to the measurement and analysis of line strengths for 
this sample of BCGs. 

This paper is organised as follows. Section 2 details the 
sample and the selection criteria, followed by a description 
of the observations and data reductions in Section 3. The 
kinematical measurements are described in Section 4 and 
the individual galaxy kinematic profiles and notes are pre- 
sented in Section 5. Section 6 discusses the kinematic prop- 



erties of BCGs as a class, compared to other Hubble types. 
Conclusions and future work are given in Section 7. 



2 SAMPLE 

The initial intention of this project was to study cD galaxies 
in particular. However, the confusing regarding the classifi- 
cation of cD galaxies and the lack of deep photometry which 
would allow to conclusively distinguish between BCGs with 
and without a halo made the sample selection extremely 
challenging. Aiming at maximising the completeness of the 
cD sample and inclusion of any potential sub-populations we 
adopted rather broad selection criteria when restricting an 
intrinsic BCG sample. We expect to add deep photometric 
information in future publications which will allow to test 
the presence of a faint halo and discriminate between cD and 
non-cD BCGs. For the purpose of this publication we refer 
to our sample of 63 galaxies as a BCG sample and present 
kinematics results for 41 observed members. Two ordinary 
elliptical galaxies were also observed (one E and one E/SO) 
and their kinematics presented here. They will be used as a 
control sample in this study, where relevant. 

The sample selection combines three methods making 
the best use of available information from literature and 
astronomical databases. The three search methods are de- 
scribed below. 

Method 1. Two well-known galaxy cluster classification 
systems that distinguish clusters containing a cD galaxy in 
the centre from other galaxy clusters, are those of Rood 
& Sastry (1971, hereafter RS) and Bautz & Morgan (1970, 
hereafter BM). The RS classification is based on the pro- 
jected distribution of the brightest 10 members and the BM 
classification is based on brightness contrast between first- 
and second-ranked galaxies (i.e. the slope of the luminosity 
function at the bright end). 

Hoessel, Gunn & Thuan (1980) derived the BM types 
for all the nearby Abell clusters of galaxies, which made it 
possible to obtain a list of all nearby Abell clusters with BM 
types I and I-II (classified to contain a cD galaxy). Struble 
& Rood (1987) compiled a catalogue of the morphological 
properties of 2712 Abell clusters, derived from visual inspec- 
tion of the KPNO photographic plates. This made it possible 
to obtain a list of Abell clusters with RS cluster classification 
type cD. The NASA/IPAC Extragalactic Database (NED0 
was used to search for the brightest galaxies close to the cen- 
tre of each cluster. If the brightest member was classified as 
either a cD or a D galaxy in NED (in the NED morphology 
or notes of previous observations) , it was included in our list 
of BCGs. The RS and/or BM classifications are frequently 
used as a method to identify possible cD galaxies (for exam- 
ple Hill & Oegerle 1992; Baier & Wipper 1995; Giacintucci 
et al. 2007). 

Method 2. The second search was carried out by perform- 
ing an all-sky search in the HyperLEDA0 database for galax- 
ies with the following properties: T-type between -3.7 and 

1 http:/ /nedwww. ipac.caltech.edu/. 

2 http://leda.univ-lyonl.fr/. 
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-4.3 (the de Vaucouleurs Third Reference Catalogue (RC3) 
classifies cD galaxies as T-type = -4); apparent B-magnitudc 
brighter than 16; distance closer than 340 Mpc; and further 
than 15 degrees from the Galactic plans[f|. It was found that 
T-type alone is not a reliable classification of galaxy types 
(in both the RC3 and HyperLeda catalogues), as various low 
luminosity galaxies are also categorised as galaxies with T- 
type = -4. This prompted an absolute magnitude cutoff of 
Mb = -20 to be applied throughout the sample selection. 
The galaxies obtained with this search were added to our 
list of BCGs if they were classified as cD galaxies in NED. 

Method 3. The third search was undertaken by choos- 
ing galaxies with published surface brightness profiles that 
had a slope consistent with cD galaxies and in the case 
of cD galaxies broke the de Vaucouleurs r* law at large 
radii. Galaxies were chosen from the series of papers by 
Schombert (1986; 1987; and 1988) and from Malumuth & 
Kirshner (1985). These are all optical photometric studies 
of the brightest cluster members of clusters. Galaxies with 
the surface brightness slopes and profiles conforming to the 
cD/D criteria set in Schombert (1987) were included in our 
list since that forms the basis of most definitions of cD galax- 
ies. 

The complete list of galaxies and their properties is 
given in Table [T] The following global criteria were applied 
to the sample: apparent B-magnitude brighter than 16; dis- 
tance closer than 340 Mpc; and the absolute magnitude cut- 
off at Mb = -20. Method 1 delivered 13 galaxies, 2 delivered 
32 galaxies and method 3 contributed 29 galaxies (indicated 
in Table [1}. With 11 of the merged list being duplicates 
present from more than one method, this provided an all- 
sky list of 63 galaxies. 

In summary: These 63 galaxies were classified as cD/D 
either in NED (in the morphological classification or in the 
notes of previous observations) and/or have profiles breaking 
the A law. With the exception of NGC4946 and NGC6047, 
which are included in the sample as control galaxies but 
known to be an E and E/SO galaxy respectively, all galaxies 
have thus been previously classified in the literature to con- 
tain a halo, albeit very inhomogeneously. Despite the confu- 
sion regarding the halo classification, all these galaxies are 
the dominant central galaxies in clusters or groups. If selec- 
tion method 1 and 3 above are assumed to be more robust 
methods of identifying a BCG with a halo, then 67 per cent 
of the BCG sample are confirmed cD galaxies. 

Ten galaxies in the BCG sample are in common with 
the Fisher et al. (1995a) sample of 13 BCGs. For six Abell 
clusters, the choice of BCG differs from that by Postman 
& Lauer (1995). These are: Abell 0189, 0194, 0262, 0548, 
4038 and 2151. In each of the first five cases, the choice of 
BCG in this sample is brighter and more dominant than the 
BCG in the Postman & Lauer (1995) sample. In the case 
of Abell 2151, a comparison of the brightness of the BCG 
choices was complicated by the fact than the BCG chosen 
by Postman & Lauer (1995), NGC6041, is actually an inter- 
acting pair of galaxies. For this cluster, we followed previous 



3 An Hq value of 75 km s 1 Mpc 1 is assumed throughout this 
work. 



authors (e.g. Aragon-Salamanca et al. 1998) who have cho- 
sen NGC6034 to be the BCG as it has a very prominent cD 
halo (Schombert 1987, 1988), and is close to an X-ray peak. 

Some of our clusters are known to have substructure 
(Rines et al. 2002; Adami et al. 2005) and, in these cases, 
there may be more than one local X-ray maximum, which 
can host a dominant galaxy central to that substructure. 
Hence, according to the definition adopted here, there can 
be more than one BCG per cluster. In particular, three of 
the selected BCGs (NGC4889, NGC4874 and NGC4839) are 
hosted by the Coma cluster, all of them showing pronounced 
envelopes (Schombert 1987; Andreon et al. 1996). The first 
two are the dominant galaxies in the cluster (it is not clear 
which one is the original central galaxy of the main cluster, 
see e.g. Adami et al. 2005), while NGC4839 is the dominant 
galaxy of a group that appears to be infalling into the main 
cluster (Neumann et al. 2001). Two galaxies (NGC6173 and 
NGC6160) from the cluster Abell 2197 are also included in 
the sample. These two galaxies are believed to be the dom- 
inant galaxies of two main groups in the process of merging 
(Abell 2197E and Abell 2197W; Schombert 1987). 

An internally consistent photometric study of this BCG 
sample will be done before investigating the properties of cD 
galaxies as a separate class. Liu et al. (2008) measured the 
extended envelopes of a sample of BCGs using Petrosian 
profiles (Patel et al. 2006) . There is a distinct signature of a 
plateau in the Petrosian profiles of cD galaxies with an ex- 
tended stellar halo which is not present for normal elliptical 
galaxies. However, they find that the deviation of the surface 
brightness profiles from a single Sersic profile does not have 
any sharp transitions. Therefore, it is very difficult to unam- 
biguously separate cD from non-cD BCGs based on surface 
brightness profiles alone (Liu et al. 2008), and confirming 
the presence or absence of cD haloes depends on choices 
of profiles fitted to the surface brightness as a function of 
radius. 

In all the conclusions of this study, it will be clearly 
stated whether it applies to BCGs or includes the two nor- 
mal elliptical galaxies. The 43 galaxies discussed further in 
this first paper are indicated in the last column of Table [1] 
(numbered 1 to 5 according to observing run), and range 
from Mb = -23.54 to Mb = -20.71 for the BCGs. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 WHT observations 

Spectroscopic observations of 10 of the 43 galaxies presented 
here (listed in Table [5J) were carried out during the period 
23 to 26 June 2006, using the 4.2m William Herschel tele- 
scope (WHT) equipped with the ISIS double spectrograph 
mounted at the //ll Cassegrain focus. The Marconi2 CCD 
was employed on the red arm and the EEV12 on the blue. 
The spatial scale of both CCDs was 0.4 arcsec/pixel (with 
2x2 binning) and the slit width was 1 arcsec. Compari- 
son spectra were provided by Cu/Ar and Cu/Ne arcs, and 
flat-field illumination by a tungsten lamp. 

Two different dichroics were used. The 5300 A dichroic 
was used for three nights and delivered an unvignetted spec- 
tral range of 3900 - 5460 A in the blue arm and 5730 - 6960 
A in the red arm. The 6100 A dichroic was used for one 
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Object a(J2000) 5(J2000) Cluster NED-type T-typc PA m B M B Ref Run 

Name degrees 



IC1565 





39 


26.3 


06 


44 


03 


A0076 


cD;E 


-3.9±1 


.2 


- 


14.41 


-21.96 


2 




UGC00579 





56 


16.2 


-01 


15 


22 


A0119 


cD;E 


-4.8±0 


4 


41 


14.40 


-22.34 


3 


6 


ESO541-013 


1 


02 


41.8 


-21 


52 


55 


A0133 


cD;E + 3 pec 


-4.0±0 


.7 


16 


14.71 


-22.54 


2 


6 


PGC004072 


1 


08 


50.8 


-15 


24 


31 


A0151 


D 


-4.2±1 


.0 


83 


14.75 


-22.37 


2 


6 


IC1633 


1 


09 


55.6 


-45 


55 


52 


A2877 


cD;E + l 


-3.9±0 


.5 


97 


12.57 


-22.60 


2, 3 


3 


IC1634 


1 


11 


02.9 


17 


39 


46 


A0154 


MLT SYS 


-4.1±1 


.0 


153 


15.27 


-22.64 


1 




IC1695 


1 


25 


07.6 


08 


41 


58 


A0193 


cD;S? 


- 




103 


14.87 


-21.98 


1 




NGC0533 


1 


25 


31.5 


01 


45 


33 


A0189 


cD;E3: 


-4.8±0 


.6 


50 


12.45 


-22.25 


3 


6 


NGC0541 


1 


25 


44.3 


-01 


22 


46 


A0194 


cD;SO-: 


-3.9±0 


.9 


69 


13.09 


-21.62 


2 


6 


IC1733 


1 


50 


42.9 


33 


04 


55 


A0260 


cD;E: 


-4.9±0 


.8 


- 


14.14 


-22.16 


1 




NGC070S 


1 


52 


46.3 


36 


09 


12 


A0262 


cD;E 


-4.9±0 


.4 


39 


13.63 


-20.98 


3 




UGC02232 


2 


46 


03.9 


36 


54 


19 


A0376 


cD;E/D 


-3.9±1 


.3 


- 


15.50 


-21.64 


2 


2 


NGC1129 


2 


54 


27.4 


41 


34 


46 


AWM7 


E;BrClG 


-4.7±1 


.2 


73 


13.34 


-21.67 


3 




UGC02450 


2 


58 


57.8 


13 


34 


59 


A0401 


cD 


2.8±5 


.0 


29 


15.24 


-23.50 


1 




NGC1275 


3 


19 


48.2 


41 


30 


42 


A0426 


cD;pec;NLRG 


-2.2±1 


.7 


110 


12.55 


-22.69 


3 




NGC1399 


3 


38 


29.0 


-35 


26 


58 


RBS454 


cD;El pec 


-4.5±0 


.5 


- 


10.42 


-20.81 


3 


3 


ESO303-005 


4 


13 


58.8 


-38 


05 


50 


RBS512 


cD? 


-4.0±0 


.5 


- 


15.38 


-21.51 


2 


3 


MCG-02-12-039 


4 


33 


37.8 


-13 


15 


40 


A0496 


cD;E + ? 


-3.9±0 


.9 


180 


13.95 


-22.50 


2, 3 


3 


ESO202-043 


4 


37 


47.6 


-51 


25 


23 


A S0479 


E+ 


-3.8±0 


.8 


133 


14.45 


-21.72 


2 


6 


ESO552-020 


4 


54 


52.3 


-18 


06 


53 


CID 28 


cD;E+ 


-3.9±0 


.7 


148 


13.54 


-22.51 


2 


3 


NGC1713 


4 


58 


54.6 


00 


29 


20 


CID 27 


cD;E + 


-4.3±0 


.6 


39 


13.88 


-20.71 


2 


2 


UGC03197 


4 


59 


55.8 


80 


10 


43 


A0505 


cD 


-4.9±0 


.9 


85 


15.25 


-22.26 


3 




ESO488-027 


5 


48 


38.5 


-25 


28 


44 


A0548 


cD;E + l 


-3.8±0 


.6 


68 


14.22 


-22.22 


2 


3 


PGC025714 


9 


08 


32.4 


-09 


37 


47 


A0754 


D 


-4.2±1 


.4 


122 


14.32 


-23.04 


2 


6 


PGC026269 


9 


18 


05.7 


-12 


05 


44 


A0780 


(R')SAO-:;BrClGSy3 


-2.6±1 


.0 


133 


14.38 


-22.84 


1 


3 


NGC2832 


9 


19 


46.9 


33 


44 


59 


A0779 


E+2:;cD 


-4.3±0 


.6 


172 


12.79 


-22.38 


2, 3 


2 


UGC05515 


10 


13 


38.3 


00 


55 


32 


A0957 


E+pec: 


-4.0±0 


.7 


83 


14.49 


-22.28 


2, 3 


5 


PGC030223 


10 


20 


26.6 


-06 


31 


35 


A0978 


D 


-3.9±1 


.2 


1 


15.40 


-21.85 


2 


6 


NGC3311 


10 


36 


42.9 


-27 


31 


37 


A1060 


cD;E + 2 


-3.3±1 


.3 


- 


12.69 


-21.33 


1 


5 


NGC3842 


11 


44 


02.2 


19 


56 


59 


A1367 


E;BrClG 


-4.8±0 


.7 


177 


12.80 


-22.18 


3 


1 


PGC044257 


12 


57 


11.4 


-17 


24 


36 


A1644 


cD;E+4 


-4.1±1 


.0 


44 


13.50 


-23.54 


1, 2 


5 


NGC4839 


12 


57 


24.2 


27 


29 


54 


A1656 


cD;SAO 


-4.0±0 


.8 


64 


13.05 


-22.26 


2, 3 


4 


NGC4874 


12 


59 


35.5 


27 


57 


36 


A1656 


cD;Di 


-3.7±0 


.9 


45 


12.75 


-22.50 


2, 3 


1 


NGC4884(9) 


13 


00 


07.9 


27 


58 


41 


A1656 


cD;E4;Db 


-4.3±0 


.5 


82 


12.48 


-22.56 


2 


1 


NGC4946* 


13 


05 


29.4 


-43 


35 


28 


A3526 


E+? 


-3.8±0 


.8 


135 


13.41 


-20.27 




5 


ESO444-046 


13 


27 


56.9 


-31 


29 


44 


A3558 


cD;E+4 


-3.7±1 


.0 


161 


14.07 


-22.86 


2, 3 


5 


LEDA094683 


13 


53 


06.4 


05 


08 


59 


A1809 


cD 


-2.3±5 


.0 


46 


15.30 


-22.35 


1 


5 


GSC555700266 


14 


01 


36.4 


-11 


07 


43 


A1837 


cD 


- 




- 


14.59 


- 


1 


5 


NGC5539 


14 


17 


37.8 


08 


10 


46 


A1890 


cD 


0.7±5 


.0 


39 


14.92 


-22.62 


1 




IC1101 


15 


10 


56.1 


05 


44 


41 


A2029 


cD;SO-: 


-2.9±1 


.1 


25 


15.10 


-23.05 


3 


5 


UGC09799 


15 


16 


44.6 


07 


01 


16 


A2052 


cD;E 


-4.7±1 


.3 


30 


14.34 


-21.84 


3 




PGC054913 


15 


23 


05.3 


08 


36 


33 


A2063 


cD;S 


-3.8±1 


.6 


- 


14.71 


-21.46 


2 




UGC09958 


15 


39 


39.1 


21 


46 


58 


A2107 


cD;SO-: 


-2.9±0 


.7 


98 


14.76 


-22.00 


3 




UGC10012 


15 


44 


59.0 


36 


06 


36 


A2124 


cD;E 


-4.9±0 


.6 


143 


14.54 


-23.15 


3 




UGC10143 


16 


02 


17.3 


15 


58 


28 


A2147 


cD;E + 


-4.0±0 


.8 


12 


14.27 


-21.97 


2, 3 


1 


NGC6034 


16 


03 


32.1 


17 


11 


55 


A2151 


E+ 


-4.0±0 


.4 


59 


14.58 


-21.60 


2, 3 


1 


NGC6047* 


16 


05 


09.1 


17 


43 


48 


A2151 


E+ 


-3.5±0 


.8 


117 


14.58 


-21.45 




1 


NGC6086 


16 


12 


35.6 


29 


29 


06 


A2162 


cD;E 


-4.8±0 


.6 


3 


13.81 


-22.19 


3 


4 


NGC6160 


16 


27 


41.1 


40 


5 5 


37 


A2197 


cD;E 


-4.8±0 


.6 


72 


14.17 


-21.70 


3 


4 


NGC6166 


16 


28 


38.5 


39 


33 


04 


A2199 


cD;E 


-4.2±1 


.4 


38 


12.88 


-22.92 


2, 3 


1 


NGC6173 


16 


29 


44.9 


40 


48 


42 


A2197 


cD;E 


-4.8±0 


.6 


138 


13.13 


-22.58 


3 


4 


NGC6269 


16 


57 


58.1 


27 


51 


16 


AWM5 


cD;E 


-4.7±0 


.7 


80 


13.38 


-22.98 


3 


4 


IC4765 


18 


47 


18.5 


-63 


19 


50 


A S0805 


cD;E+4 


-3.9±0 


.5 


123 


12.46 


-21.97 


2 


5 


NGC7012 


21 


06 


45.8 


-44 


48 


49 


A S0921 


E+4 pec 


-3.9±0 


.5 


100 


13.93 


-21.79 


2 


3 


ESO146-028 


22 


28 


51.1 


-60 


52 


55 


RXCJ2228.8-6053 


E+3 


-3.8±0 


.8 


154 


13.94 


-22.54 


2 


5 


ESO346-003 


22 


49 


22.0 


-37 


28 


20 


A S1065 


E + 2 pec: 


-3.8±0 


.8 


118 


14.00 


-21.59 


2 


5 


NGC7597 


23 


18 


30.3 


18 


41 


20 


A2572 


cD;S? 


-2.0±0 


.8 


133 


15.00 


-21.52 


1 


2 


NGC7647 


23 


23 


57.4 


16 


46 


38 


A2589 


cD;E 


-4.7±0 


.8 


174 


14.61 


-21.93 


3 


1 


NGC7649 


23 


24 


20.1 


14 


38 


49 


A2593 


cD;E 


-4.6±1 


.1 


78 


15.65 


-20.99 


3 


2 


PGC071807 


23 


35 


01.5 


27 


22 


20 


A2622 


cD 






138 


15.28 


-22.05 


1 


4 


NGC7720 


23 


38 


29.5 


27 


01 


51 


A2634 


cD;E + pcc: 


-4.3±0 


.6 




13.43 


-22.63 


2 


1 


IC5358 


23 


47 


45.0 


-28 


08 


26 


A4038 


cD;E+4 pec 


-3.8±1 


.1 


114 


13.64 


-22.00 


2 


3 


NGC7768 


23 


50 


58.3 


27 


08 


51 


A2666 


cD;E 


-4.9±0 


.5 


55 


13.28 


-22.39 


3 


1 


PGC072804 


23 


54 


13.7 


-10 


25 


09 


A2670 


cD;XBONG 


-4.4±2 


.4 


76 


15.32 


-22.82 


1 


5 


ESO349-010 


23 


57 


00.7 


-34 


45 


33 


A4059 


cD;E+4 


-3.7±1 


.2 


155 


14.18 


-22.69 


2 


3 



Table 1. The new sample of BCGs, including the two ordinary ellipticals (marked *). Columns 2, 3, 6 — 9 arc from the HyperLeda 
catalogue, and 4 and 5 from NED. PA is the position angle of the galaxy major axis. Ref 1,2 and 3 correspond to the three search 
methods for objects (column 10). Run (column 11) 1=WHT, 2=Gemini N 2006B, 3=Gcmini S 2006B, 4=Gemini N 2007A, 5=Gemini S 
2007A, 6=Gcmini S 2007B (The data from Run 6 are not presented in this paper). Besides the galaxies classified as cD/D in the NED 
morphology (column 5); NGC6034, UGC05515 and NGC1129 were classified as cD/D in Schombert (1987); NGC7012, ESO146-028, 
ESO346-003, ESO202-043, NGC3842, PGC026269 and IC1634 were classified as cD/D in the NED notes by previous observations. 
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night, delivering an unvignetted spectral range of 4500 - 
6120 A in the blue arm and 7950 - 9600 A in the red arm. 
The exposure times are listed in Table [2] This combination 
made it possible to observe 29 absorption indices consisting 
of the original 21 Lick indices from CN to Ti02 (Burstein et 
al. 1984; Faber et al. 1985; Worthey et al. 1994), the four H<5 
and H7 indices contributed by Worthey & Ottaviani (1997), 
the near-IR Call triplet at ~ 8600 A and the Mgl index. 
This collection of 29 indices will be referred to as the Lick 
indices in this study. 

Dispersion in the blue arm was 0.90 A/pixel (grating 
R600B) and 0.88 A/pixel (grating R600R) in the red. The 
spectral resolution ranged between 3.3 and 4.0 A, depend- 
ing on the dichroic and arm used. Seeing was typically better 
than 1 arcsecond, and all necessary calibration frames were 
obtained. In most cases, the slit was placed on the major 
axis of the galaxy. Exceptions occurred if there was another 
early-type galaxy (unrelated to the sample) in close proxim- 
ity, in which case the slit was positioned to go through the 
centres of both galaxies, to allow a control sample of ellipti- 
cals with the same observational setup (two control sample 
ellipticals were observed but are not used in this present 
paper). Five spectrophotometric standard stars were also 
observed for flux calibration. 

A total of 22 Lick calibration stars (G and K type stars 
which are considered the main contributors to the visible 
light in early-type galaxies) were observed with the 5300 A 
dichroic and 10 with the 6100 A dichroic. A neutral density 
filter was used for the Lick stars with both dichroics. These 
stars are used as templates for velocity dispersion measure- 
ments as well as to transform the line-strength indices to 
the Lick system. Before and/or after every standard star or 
object spectrum, an arc spectrum and flat-field frame were 
obtained in the red and blue arm at each telescope pointing, 
to allow accurate wavelength calibration as well as removal 
of the response function of the dichroics. 

3.2 Gemini observations 

The Gemini Multi-object Spectrograph (GMOS) observa- 
tions were performed in queue scheduling mode. The B600 
grating was used (2.7 A resolution) with a 0.5 arcsec slit 
width to match the observations of the Lick stars described 
below. Dispersion was 0.914 A/pixel, and the spatial scale 
was 0.146 arcsec/pixel (with 2x2 binning). Ten galaxies 
were observed at Gemini South in semester 2006B (August 
2006 to January 2007) and five galaxies at Gemini North in 
the same observing semester. A further 12 galaxies were ob- 
served at Gemini South and six at Gemini North in 2007A 
(February 2007 to August 2007). 

Bryan Miller et al. (private communication) performed 
a series of observing programmes from 2002 to 2004 in which 
they acquired long-slit spectroscopic data of "Lick Stars". 
These are well suited as calibration data for our observ- 
ing programmes, and our instrumental set-up corresponded 
to the set-up used for observations of the Lick stars. Spec- 
tral dithering was carried out with two central wavelengths 
at 5080 A and at 5120 A for the Lick stars as well as the 
galaxy data described here. This is done in order to obtain 
the full spectrum uninterrupted by the gaps between the 
three CCDs. The three CCDs mosaiced together delivered a 
spectral range from 3700 to 6500 A. Exposure times were 2 



x 900 seconds at each position, which amounts to an hour 
(mostly along the major axis) for each galaxy, except in the 
case of ESO444-046 where only one 900 second exposure at 
each position was taken. Imaging sequences (3 x 20 seconds 
with the g and r filter) were performed to allow insight into 
the possible presence and location of dust lanes or gas rel- 
ative to the slit position. Standard stars were observed for 
relative flux calibration. The observations were mostly exe- 
cuted in dark time. The seeing typically ranged from 0.6 to 

1.2 arcsec, and very rarely went up to 2.2 arcsec. Calibra- 
tion arc (CuAr) and fiat field spectra were also frequently 
observed at the two different central wavelengths. General 
bias frames and imaging flat field frames were also obtained. 

For the Gemini data, information on the slit placements 
during observations was used to reconstruct the slit position 
on images of the galaxies. The slit location on the targets 
was determined using the slit acquisition images and the tar- 
get positions from the imaging sequences. The images were 
taken with the same telescope pointing, and used to measure 
the offsets of the slit position from the galaxy centre (Table 
In the few cases where the same galaxy was observed on 
two different nights, this was done for both nights. Offsets 
from the galaxy centres are shown in Table [3] The major- 
ity of the galaxies were centred correctly or with very small 
offsets. There are no offsets indicated for the WHT obser- 
vations, because these were carried out in visiting observer 
mode, and all offsets were zero. 

3.3 Data reduction 

All the WHT data reductions were performed with the 
iRAiJ^and STARLlNpJf] reduction packages. The blue and red 
spectra and spectra taken using the 5300 and 6100 A 
dichroics were treated in similar parallel streams of data 
reduction. Overscan correction, ADU conversion, bias sub- 
traction and trimming of the image edges were done first. 
Flat field spectra were normalised by correcting for the lamp 
response function and all star and galaxy spectra were ap- 
propriately flat-fielded. The dichroics introduced an inter- 
mediate frequency pattern which varied with the position of 
the telescope. This was successfully removed by flat-fielding, 
using the flat-field images taken directly before and/or af- 
ter the science frames (i.e. at the same telescope pointing). 
It was necessary to make an illumination correction using 
normalised twilight flat field spectra to correct for spatial 
gradients in the star and galaxy spectra. Full 2D wavelength 
calibration was performed and star and galaxy spectra were 
transformed to equal linear spectral bin sizes. 

All the flux calibration and Lick stars were sky- 
subtracted and S-distortion-corrected before ID spectra 
were extracted. All ID Lick star spectra were flux calibrated 
and corrected for atmospheric extinction. 

A cosmic ray cleaning algorithm was applied to all the 
2D galaxy spectra before they were background subtracted 

4 I RAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy Inc., under cooperative agreement with 
the National Science Foundation. 

STARLINK is an interactive reduction and analysis facility that 
was developed in the UK, funded by the PPARC/STFC Research 
Council. 
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Object Exposure Exposure Slit Major Axis r e e a c Fraction 

(5300 A) (6100 A) PA MA times 

seconds seconds degrees degrees arcsec arcsec a e 



NGC3842 


6x1200 


2x900 


5 


177 


16.4 


0.29 


16.3 


1.22 


NGC4874 


6x1200 


3x900 


82 


45 


20.0 


0.00 


20.0 


0.60 


NGC4889 


6x1200 


2x900 


80 


82 


25.7 


0.34 


25.7 


0.89 


NGC6034 


6x1200 


3x900 


50 


59 


8.2 


0.27 


8.1 


1.23 


NGC6047 


3x1200 




88 


117 


7.4 


0.27 


7.0 


0.61 


NGC6166 


6x1200 


3x900 


35 


38 


22.6 


0.26 


22.6 


0.22 


NGC7647 


6x1200 


3x900 


174 


174 


11.2 


0.49 


11.2 


1.21 


NGC7720 


6x1200 


3x900 


10 




9.3 


0.19 


9.3 


0.65 


NGC7768 


3x1200 


3x900 


60 


55 


14.7 


0.19 


14.7 


0.82 


UGC10143 


6x1200 


3x900 


12 


12 


10.5 


0.41 


10.5 


0.65 



Table 2. Exposure times in each of the two dichroics and position angles for the WHT BCGs. The PA is given as deg E of N. The 
half-light radii (r e ) were calculated from the 2MASS catalogue, except NGC6166 and NGC7720 for which they were calculated with the 
VAUCOUL task in the reduceme package. The last column lists the fraction of the effective half-light radii spanned by the radial profiles 
measured in this work. Table [T] contains more detailed properties for the galaxies. 



Object Slit Major Axis r e e a e Fraction Offset 

PA MA arcsec arcsec a e arcsec 







2006B South 


















ESO303-005 


55 




9.1 





.25 


9. 


1 





.50 





.70 


ESO349-010 


14 


155 


15.3 


0. 


.52 


12. 


3 


1 


.18 





.30 


ESO488-027 


88 


68 


10.3 


0. 


.15 


10. 


1 





.94 


2. 


.20 


ESO552-020 


148 


148 


18.3 





.43 


18. 


3 


0. 


.76 


0. 


.64 


IC1633 


97 


97 


23.9 





.17 


23. 


9 


0. 


.75 


0. 


.00 


IC5358 


40 


114 


17.4 


0. 


.60 


8. 


3 


1. 


.11 





.26 


MCG-02-12-039 


166 


180 


15.2 


0. 


.19 


15. 


1 


1 


.13 





.64 


NGC1399 


222 




42.2 





.06 


42. 


2 


0. 


.41 


0. 


.45 


NGC7012 


289 


100 


15.8 





.44 


15. 


6 


0. 


.96 


0. 


.00 


PGC026269 


313 


133 


7.3 


0. 


.00 


7. 


3 





.89 





.00 






2006B North 


















NGC1713 


330 


39 


15.5 





.14 


14. 





1. 


.07 


0. 


.00 


NGC2832 


226 


172 


21.2 





.17 


19. 


6 


1. 


12 


0. 


.50 


NGC7597 


46 


133 


9.4 


0. 


.00 


9. 


A 





.64 





.00 


NGC7649 


78 


78 


12.8 


0. 


.31 


12. 


.8 


0. 


.82 





.20 


UGC02232 


60 




9.7 





.00 


9. 


7 


0. 


.82 


0. 


.00 






2007A South 


















ESO146-028 


154 


154 


12.4 


0. 


.35 


12. 


.4 


0. 


.48 





.15 


ESO346-003 


297 


118 


11.9 


0. 


.18 


11. 


9 





.29 





.15 


ESO444-046 


196 


161 


16.0 





.29 


14. 


9 


0. 


.24 


0. 


,15 


GSC555700266 


204 




10.6 





.30 


10. 


6 


0. 


.47 


0. 


.00 


IC1101 


204 


25 


12.2 


0. 


.50 


12. 


2 


0. 


.39 





.20 


IC4765 


287 


123 


28.0 


0. 


.46 


27. 


1 





.13 


0. 


.75 


LEDA094683 


226 


46 


7.3 





.33 


7. 


3 


0. 


.54 


2. 


00 


NGC3311 


63 




26.6 





.17 


26. 


6 


0. 


.53 


0. 


.40 


NGC4946 


134 


135 


17.0 


0. 


.20 


17. 








.74 





.00 


PGC044257 


267 


41 


15.3 


0. 


.63 


10. 


5 


0. 


.38 


0. 


.00 


PGC072804 


76 


76 


7.6 





.16 


7. 


6 


0. 


.52 





.00 


UGC05515 


293 


83 


12.0 





.13 


11. 


.8 


0. 


.38 


0. 


.00 






2007A North 


















NGC4839 


63 


64 


17.2 


0. 


.52 


17. 


2 


0. 


.61 





.00 


NGC6086 


270 


3 


11.2 





.29 


8. 


1 


0. 


.51 


0. 


.10 


NGC6160 


140 


72 


13.9 





.47 


9. 





0. 


.41 


0. 


.00 


NGC6173 


139 


138 


15.0 


0. 


.26 


15. 








.47 


0. 


.10 


NGC6269 


306 


80 


14.1 


0. 


.20 


13. 


1 





.50 





.35 


PGC071807 


138 


138 


6.7 





.70 


6 


7 


0. 


.71 


0. 


.10 



Table 3. Galaxies observed with the Gemini North and South telescopes. Table [T] contains more detailed properties for the galaxies. 
The PA is given as deg E of N. The half-light radii (r e ) were calculated from the 2MASS catalogue, except ESO303-005 for which it 
was calculated with the VAUCOUL task in the REDUCEME package. The second last column lists the fraction of the effective half-light radii 
spanned by the radial profiles measured in this work, and the last column lists the offsets as described in the text. 
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and S-distortion corrected. For each individual galaxy frame, 
a sky region consisting of 40 rows was chosen on each side 
of the galaxy and averaged together to obtain a ID sky 
spectrum. The sky spectrum was grown into a 2D spectrum 
and subtracted from the data frame. Flux calibration and 
extinction corrections were also applied to the 2D spectra. 
The galaxy frames were then co-added (also removing any 
remaining cosmic ray effects in the process) to form a single 
frame for each galaxy in each of the blue and red arms and 
for each of the two dichroics. 

The basic Gemini data reductions were performed with 
the Gemini-specific GMOS data reduction package that can 
be added onto the standard iraf reduction package. Spectra 
with different central wavelengths (5080 and 5120 A) and 
spectra taken with GMOS-North and GMOS-South were 
treated in similar parallel streams of data reduction. The 
raw Lick star data were obtained from the Gemini Science 
Archiv^f). The basic data reductions were conducted in a 
similar manner as for the WHT data. For each of the two 
central wavelength settings, a bad pixel spectrum was made 
in order to avoid interpolation between the gaps. When the 
spectra of the two different wavelength settings were com- 
bined, these bad pixel masks were used to ensure that only 
true data were averaged. 

Proctor et al. (2005) and Brough et al. (2007) also used 
GMOS data to measure Lick indices, and found negligible 
scattered light that influenced their results. Nevertheless, as 
a test of the effect of scattered light in the GMOS spectro- 
graph on our data, bright skylines were used to extract flat 
field profiles that could be compared to the shape of the 
dome flats that were used. It was found that scattered light 
effects did not influence the regions where galaxy data were 
extracted or the regions used for background subtraction. 

An error analysis was completed to account for both 
Poisson (photonic) noise and systematic errors in the data. 
Variance frames were created from the data frames after 
bias subtraction and ADU conversion, and reduced in a way 
similar to the data frames, but the arithmetic manipulations 
were performed by incorporating the usual error propagation 
laws. Thus, the variance frames represent the square of the 
Poisson noise at the start of this reduction process, and all 
other systematic noise effects are added onto the variance 
frames in quadrature. 

Flat-fielding produced negligible errors with the excep- 
tion of the WHT flat-field spectra taken in the upper wave- 
length region (red spectra of dichroic 6100 A on night four) 
which contained fringe patterns. This contributed about four 
per cent to the errors and these were appropriately propa- 
gated in the variance frame. Errors contributed by the illu- 
mination correction proved negligible. Wavelength calibra- 
tion was accurate to within 0.1 A in all cases. After that, 
the variance frames were also S-distortion corrected, follow- 
ing the procedure used for the corresponding data frames. 
The sky subtraction process added to the uncertainty in the 
data and was properly taken into account by calculating the 
following for each variance frame: 

[ASky(\)] 2 = [JL AA(X)} 2 , (1) 



http: / /www4.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/gsa/ 



where ASky is the error contribution of the sky subtraction 
being taken into account, N is the number of spatial cross 
sections averaged for the sky determination and A is the 
Poisson error on the average one line sky spectrum. This 
ID spectrum was then grown into a 2D spectrum with the 
proper dimensions and added onto the variance frames in all 
cases. 

Relative flux calibration was accurate to within five to 
10 per cent, and the variance frames were also co-added 
so that they could be treated alongside the galaxy data in 
further analysis. For the Gemini data, variance frames were 
co-added for each galaxy by dealing with the gaps in the 
same manner as the data (described above) to avoid any 
interpolation. The fact that only N /2 frames were added in 
the regions where the gaps occur, instead of N individual 
frames per galaxy in all the other regions, was taken into 
account in the variance frame combination procedure. At 
the end of this reduction process, each 2D galaxy data frame 
had an associated 2D variance array. 

As mentioned, some of the galaxies were observed at a 
position angle (PA) other than the major axis (MA) , due to 
the two control elliptical galaxies in the WHT observations, 
and due to the availability of suitably bright guide stars 
in the Gemini observations. This had an effect on the rela- 
tive half-light radii (r c ) which had to be taken into account. 
Tables [2] and [3] list the r c along the semimajor axis of the 
galaxies, computed from the 2MASS A"-band 20th magni- 
tude arcsec -2 isophotal radius using the formula by Jarrett 
et al. (2003): logr e ~ logrK20 — 0.4. Using the ellipticities 
(e) of the galaxies (data from NED), the effective half-light 
radii at PA were computed according to 



c 1-e | cos(\ PA- MA |) |' y 1 

The fraction of a a which the radial profiles measured in this 
work spans is also listed in Tables [2] and [3] For those galaxies 
for which a MA angle is not available in HyperLeda, the PA 
is assumed to be the same as the MA (as intended during 
the observation preparations). 



4 CENTRAL KINEMATIC AL 
MEASUREMENTS 

The central radial velocity and velocity dispersion values of 
the galaxies were derived for regions with the size of a c /8, 
and with the galaxy centres defined as the luminosity peaks. 
Thus, the central values are for an aperture of 1 x a c /8 
arcsec 2 for the WHT data and 0.5 x o e /8 arcsec 2 for the 
Gemini data. 

The Lick stars were shifted to zero radial velocity be- 
fore using them as templates. Radial velocities and velocity 
dispersions of the galaxies were measured using the movel 
and OPTEMA algorithms described by Gonzalez (1993) and 
integrated into the REDmE data reduction packagq^ This 
routine is an improvement of the classic Fourier Quotient 



REDmE is an astronomical data reduction package, specially 
devoted to the analysis of long-slit spectroscopy data. It was 
developed by N. Cardiel and J. Gorgas (Cardiel 1999). 
http://www.ucm.es/info/Astrof/software/reduceme/reduceme.html 
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Figure 1. A typical example of the calculation of the derived 
kinematics for NGC7768 with the MOVEL algorithm. In the lower 
panel, the smooth line corresponds to the optimal template spec- 
trum, plotted together with the galaxy spectrum which has been 
shifted and broadened to the derived parameters. The upper panel 
shows the residuals of fit (in relative flux units), with the dotted 
lines indicating the positions of typical emission lines. This illus- 
tration spans 4900 to 5400 A. 
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Figure 2. Velocity dispersion results compared to literature val- 
ues. The latter values are the most recent ones from HyperLeda 
as described in the text. 



Technique described by Sargent et al. (1977). It is an it- 
erative procedure in which a galaxy model is processed in 
parallel with the galaxy spectrum. In this way, a compari- 
son between the input and recovered broadening functions 
for the model allows one to correct the galaxy power spec- 
trum for any imperfections of the data handling procedures 
in Fourier Space. The procedure uses an algorithm that ob- 
tains the linear combination of individual templates that 
best matches the observed galaxy spectrum. The best tem- 
plate is scaled, shifted and broadened to match the galaxy 



kinematic parameters and is feeded to the movel algorithm. 
The output kinematic parameters are then used to create 
and improve composite templates and the process continues 
until it converges. This OPTEMA algorithm is able to over- 
come the typical template mismatch problem (see van der 
Marel & Franx 1993) by constructing an optimal template 
as a linear combination of stellar spectra of different spectral 
types and luminosity classes for each galaxy (Gonzalez 1993; 
Sanchez-Blazquez et al. 2006). A new optimal template was 
constructed with each spatial bin in the kinematical mea- 
surements. 

Skyline residuals and emission lines were masked. Emis- 
sion lines were detected in 13 galaxies: ESO349-010, MCG- 
02-12-039, NGC0541, NGC1713, NGC3311, NGC4874, 
NGC4946, NGC6166, NGC6173, NGC7012, NGC7649, 
NGC7720 and PGC044257. The detection thresholds and 
emission-line correction will be described in a future paper 
devoted to the stellar populations. 

For the Gemini data, 20 stars were used to make the 
optimal template (excluding the stars not of G or K spectral 
type) , and 22 stars for the WHT data. Figure [1] shows a 
typical fit between the observed central spectrum of a galaxy 
and the optimal template for that galaxy corrected by the 
derived kinematic parameters |f] 

The errors were computed through Monte Carlo nu- 
merical simulations and using the appropriate error spectra. 
In each simulation, gaussian noise was added to the galaxy 
spectra in accordance with the error spectra and the radial 
velocity and velocity dispersion derived. Errors in the radial 
velocities and velocity dispersions were calculated as the un- 
biased standard deviation of the different solutions. The er- 
rors computed in this manner are expected to incorporate all 
the uncertainties, from the basic reduction process (which is 
accounted for in the error spectra) to the final measurement 
of the parameters. 

When comparing measured central velocity dispersion 
values with those from the literature, emphasis must be 
placed on the fact that published values are measurements 
from within different apertures, which can have a systematic 
effect on the values. Accurate comparison is further limited 
because of different effective spatial resolution, data of dif- 
ferent quality and different analysis techniques. 

Figure [2] shows the rough comparison (not taking into 
account the above differences) between the velocity disper- 
sion values from HyperLeda (if available) and those derived 
in this study. For most of the galaxies, more than one ve- 
locity dispersion values are listed in HyperLeda (in which 
case the most recent was used for this comparison). In the 
case of NGC1399, for example, published values range be- 
tween 233 km s~ and 420 km s , which shows the effect 
of different apertures and measurements have (this can be 
seen in the steep change in the measured velocity dispersion 
profile of NGC1399 in Figure [5|. The absolute mean differ- 
ence between the values measured here and the literature is 
|<r(here) — a(lit, uncorrected)) = 33 ± 29 km s -1 , which is 
acceptable for the purposes of this paper. 



8 The derived recession velocities were automatically corrected 
for the standard special relativistic correction: z = (1 + ^ B ^)"f — 1 



where 7 : 



Vc 2 -V" 2 



, and V nr is the non-relativistic radial velocity 



and c the speed of light. 
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Galaxy a V Vmax log( Vmax / cr ) * 





Observed 


Observed 


km s 






ESO146-028 


299±3 


12061±3 


14±18 


-1 


,2±0. 


ESO303-005 


276±5 


14526±4 


23±17 


-0 


8±0. 


ESO346-003 


226±4 


8459±9 


51±13 


-0 


3±0 


ESO349-010 


282±3 


14420±3 


46±19 


-0 


.8±0. 


ESO444-046 


292 ±3 


13741±3 


26±16 


-0 


,9±0. 


ESO488-027 


248 ±2 


11754±4 


53±12 


-0 


,3±0 


ESO552-020 


229±3 


9128±5 


20±19 


-1 


,0±0 


GSC555700266 


312±9 


20059±6 


49±14 


-0 


,6±0. 


IC1101 


378±5 


22585±9 


41±31 


-1 


,0±0 


IC1633 


400±2 


7061±2 


20±17 


-1 


,0±0 


IC4765 


286±5 


4403±6 


50±28 


-0 


,8±0. 


IC5358 


243±3 


8539±9 


55±13 


-0 


,7±0 


Lcda094683 


332±5 


22748±5 


32±26 


-0 


9±0. 


MCG-02-12-039 


271±5 


9648±3 


59±28 


-0 


,3±0 


NGC1399 


371±3 


1406±2 


24±11 


-0 


,6±0 


NGC1713 


251±2 


4472 ±2 


30±11 


-0 


,5±0 


NGC2832 


364±4 


6841±4 


64±22 


-0 


,4±0 


NGC3311 


196±2 


3709±2 


25±9 


-0 


.6±0. 


NGC3842 


287±5 


62H±2 


11±13 


-1 


,2±0 


NGC4839 


278±2 


7380±2 


44±14 


-0 


,8±0 


NGC4874 


267±4 


6953±2 


15±12 


-0 


,4±0 


NGC4889 


380±4 


6266±5 


40±19 


-0 


,8±0 


NGC4946 


201±2 


3148±9 


62±9 


-0 


,2±0 


NGC6034 


325±4 


10113±28 


134±15 


-0 


,2±0 


NGC6047 


230±2 


9031±8 


59±16 


-0 


,4±0 


NGC6086 


318±5 


9483±2 


13±18 


-1 


,2±0 


NGC6160 


266±3 


9428±3 


12±18 


-1 


,3±1 


NGC6166 


310±10 


9100±2 


31±20 


-0 


.7±0 


NGC6173 


304±3 


8766±2 


24±16 


-0 


.9±0. 


NGC6269 


343±5 


10363±3 


26±35 


-0 


.8±0 


NGC7012 


240±3 


8653±3 


24±16 


-0 


.9±0 


NGC7597 


264±8 


10907±5 


63±20 





2±0 


NGC7647 


271±5 


11842±7 


36±14 


-0 


.9±0 


NGC7649 


250±5 


12258±5 


41±14 


-0 


.6±0 


NGC7720 


409±5 


9034±5 


50±18 


-0 


.6±0 


NGC7768 


272±5 


7875±24 


114±11 


-0 


.1±0 


PGC026269 


222±11 


15816±8 


51±20 





2±0 


PGC044257 


247±9 


14014±6 


20±16 


-1 


.2±1 


PGC071807 


315±3 


17911±5 


16±13 


-1 


,5±1 


PGC072804 


311±5 


22433±4 


50±24 


-0 


.4±0 


UGC02232 


314±4 


14208±3 


42±18 





0±0 


UGC05515 


362 ±4 


13103±7 


56±18 


-0 


,4±0 


UGC10143 


262±2 


10344±5 


19±11 


-1 


.1±0 



Table 4. Central velocity dispersions (cr) and radial velocities 
(V) of the galaxies. The measurement and determination of Vrn ax 
and (Vmax/ffo)* are described in Section [6] 

5 KINEMATICS OF INDIVIDUAL GALAXIES 

The galaxy spectra were binned in the spatial direction to 
ensure a minimum S/N of 20 per A in the H/3 region of 
the spectrum for kinematical measurements as a function 
of radius. Thus, the spatial cross-sections are broader with 
increasing radius from the centre of the galaxy. The radial 
kinematics (radial velocities and velocity dispersions) were 
measured with the same movel and OPTEMA algorithms in 
the REDmE package and with the same optimal Lick star 
templates as the central kinematics. The values of the kine- 
matic parameters are evaluated at the luminosity-weighted 
centres of the spatial bins used to derive the parameters. 

Dark time at both Gemini telescopes is highly oversub- 
scribed. To speed up the completion of this multi-semester 
project, we relaxed the observing conditions to accept grey 
time as well. However, as a consequence it was not possible 
to extract radial profiles to very large radii for those targets 
observed during grey time, but the central values are still ac- 
curate. Eight galaxies (three WHT and five Gemini) could 
be extracted to radii > la e and 12 galaxies (one WHT and 
11 Gemini) could only be extracted to radii < 0.5a o . The 



rotational velocity Vmax was measured as half the difference 
between the peaks of the rotation curve. 

5.1 Notes on individual objects 

This section describes the individual objects and their ra- 
dial kinematic profiles as shown in Figures [21 HI [SI El an d El 
We refer to results published prior to this study, as well as 
the results obtained here. The radial kinematic profiles of 
NGC4839, NGC4889, NGC6166 and IC1101 are compared 
to previous measurements from the literature in Appendix 

m 

ESO146-028: The radial velocity of this galaxy is mostly 
flat and the velocity dispersion profile is decreasing out- 
wards. 

ESO303-005: The radial kinematic profiles of this 
galaxy show no significant rotation or substructure. 

ESO346-003: The radial velocity profile shows clear ro- 
tation of the order of 51 km s _1 . 

ESO349-010: HST imaging of this galaxy was published 
by Laine et al. (2003). The radial velocity profile of this 
galaxy shows no significant rotation and the velocity disper- 
sion profile has a positive gradient. 

ESO444-046: No clear rotation or velocity substructure 
is visible in the kinematic profiles. 

ESO488-027: The kinematic profiles show a KDC in the 
centre of this galaxy. 

ESO552-020: The radial velocity profile is mostly flat 
and the velocity dispersion profile shows a small dip south- 
west of the centre of the galaxy which might indicate a colder 
component. 

GSC555700266: Slow rotation but no clear substructure 
can be seen in the kinematic profiles. 

IC1101: This galaxy is one of the few known BCGs that 
were previously found to have a rising velocity dispersion 
profile (Fisher et al. 1995a). This result cannot be confirmed 
here since the derived velocity dispersion profile is not mea- 
sured to the same radius as in the above mentioned study. 
No significant rotation is detected in the radial velocity pro- 
file. 

IC1633: The surface brightness profile was published by 
Schombert (1986), and HST imaging by Laine et al. (2003). 
The radial velocity profile reveals some evidence for sub- 
structure at the centre of this galaxy. 

IC4765: The radial velocity profile shows rotation in the 
centre of the galaxy which is possibly rotational substruc- 
ture. 

IC5358: Both the radial kinematic profiles show the 
clear presence of a KDC in this galaxy. 

LEDA094683: No significant rotation but some evidence 
of substructure is visible from the radial velocity profile. 

MCG-02-12-039: This BCG is a weak line-emitter. 
The surface brightness profile was published by Schombert 
(1986), and HST imaging by Laine et al. (2003). A KDC is 
detected in the radial kinematic profiles. 

NGC1399: The surface brightness profile was published 
by Schombert (1986). The velocity dispersion in the central 
region of this galaxy was found to be 388 km s _1 decreasing 
to 200 km s _1 at 10 arcseconds from the centre by Longo 
et al. (1994). The central radial velocity dispersion of this 
galaxy was measured in this study as 371 ± 3 km s _1 . Thus, 
the central velocity dispersion measured in this study cor- 
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Figure 3. Radial profiles of velocities (V) and velocity dispersion (cr). The vertical lines indicate the centre of the galaxy, and the 
horizontal line the radial velocity of the central bin. 



Loubser et al. 




r (arcsec) 

(a) IC1633. 



4460 

4440 

S" 4420 
E 

— 4400 
> 

4380 

4360 

4340 

4320 
380 

360 



r (arcsec) 

(b) IC4765. 



} 


L rffM 


} 




1 


1 kpc 



RSBO 
8F,fiO 
3540 
_ 8520 
1 8500 
> 8480 
8460 
8440 
8420 
8400 




r (arcsec) 

(c) IC5358. 





r (arcsec) 

(d) Lcda094683. 



r (arcsec) 

(c) MCG-02-12-039. 




r (arcsec) 

(f) NGC1399. 




Figure 4. Radial profiles of velocities (V) and velocity dispersion (a). The dashed lines are as in Figure[3] 
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(g) NGC6047. (h) NGC6086. (i) NGC6160. 

Figure 5. Radial profiles of velocities (V) and velocity dispersion (a). The dashed lines are as in Figure [3] 
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(g) NGC7649. (h) NGC7720. (i) NGC7768. 

Figure 6. Radial profiles of velocities (V) and velocity dispersion (cr). The dashed lines are as in Figure[3] 
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Figure 7. Radial profiles of velocities (V) and velocity dispersion (a). The dashed lines are as in Figure [3] 



16 Loubser et al. 



responds to that measured by Longo et al. (1994), and the 
steep decreasing velocity dispersion profile is confirmed as 
well. This galaxy has a flat radial velocity profile. 

NGC1713: The radial velocity profile of this galaxy re- 
veals a definite KDC in the centre of this galaxy. 

NGC2832: This galaxy is thought to be tidally interact- 
ing with its companion, NGC2831. The surface brightness 
profile were published by Schombert (1986) and Jordan et 
al. (2004), and HST imaging by Laine et al. (2003). The 
radial velocity profile shows the presence of a KDC in the 
centre of this galaxy. 

NGC3311: HST imaging was published by Laine et al. 
(2003). No clear rotation or velocity substructure can be 
seen in the radial kinematic profiles, but the velocity disper- 
sion profile has a positive gradient. 

NGC3842: This galaxy belongs to a cluster which forms 
part of the Coma supercluster. The surface brightness pro- 
file was published by Schombert (1986), and HST imaging 
by Laine et al. (2003). The radial velocity profile of this 
galaxy reveals the presence of a kinematically decoupled core 
(KDC) in the centre. 

NGC4839: The location of this cD galaxy far from the 
centre of the Coma cluster is unusual. Surface brightness 
profiles were published by Schombert (1986), Oemler (1976) 
as well as Jordan et al. (2004), who all confirm the presence 
of a very prominent cD envelope. Rotation of the order 44 
km s _1 and a KDC in the centre of the galaxy are detected. 

NGC4874: This galaxy has a large, extended envelope. 
It is the second brightest galaxy of the famous pair of cDs 
at the centre of the Coma cluster. The surface brightness 
profile was published by Peletier et al. (1990). No significant 
rotation or velocity substructure is detected for this galaxy. 

NGC4889: A large cD galaxy with a very extended en- 
velope. This is the brightest galaxy of the Coma cluster. 
The surface brightness profile was published by Peletier et 
al. (1990), and HST imaging by Laine et al. (2003). The 
radial velocity profile clearly shows a KDC. 

NGC4946: This is an ordinary elliptical galaxy. The ra- 
dial velocity profile of this galaxy reveals rotation of the 
order of 62 km s _1 along the major axis. 

NGC6034: The surface brightness profile was published 
by Schombert (1986). The radial velocity profile of this 
galaxy shows significant rotation of the order of 134 ± 15 
km s _1 along the major axis. 

NGC6047: The surface brightness profile was published 
by Schombert (1986), and Schombert (1987) classified it as 
an E/SO galaxy. The radial velocity profile of this galaxy 
shows rotation of 59 km s _1 (with the slit placement 29 
degrees from the major axis). 

NGC6086: The surface brightness profile was published 
by Schombert (1986), and HST imaging by Laine et al. 
(2003). Carter et al. (1999) found this galaxy to have a KDC. 
Even though some evidence of substructure can be seen, this 
core cannot be confirmed since the profiles measured here 
do not cover the same radial extend as those in Carter et al. 
(1999). 

NGC6160: The radial velocity profile of this galaxy re- 
veals no significant rotation or substructure. 

NGC6166: This is a classic multiple nucleus cD galaxy 
in a rich cluster. The surface brightness profile was pub- 
lished by Schombert (1986). Carter et al. (1999) found the 
velocity dispersion of this galaxy to increase from 325 km 



s _1 at the centre to 450 km s _1 at 35 arcsec along the major 
axis. They also found the galaxy to show modest major-axis 
rotation (45 km s _1 at 40 arcsec). The profiles derived here 
do not extend out to the radii necessary to confirm the pos- 
itive velocity dispersion slope found by Carter et al. (1999). 
Rotation of the order of 31 km s~ x is found in the centre of 
the galaxy although this is not a very clear rotation curve. 

NGC6173: The surface brightness profile was published 
by Schombert (1986), and HST imaging by Laine et al. 
(2003). The derived radial kinematic profiles of this galaxy 
show no significant rotation or substructure. 

NGC6269: Surface brightness profiles were published by 
Schombert (1986) and by Malumuth & Kirshner (1985). The 
radial velocity profile shows some velocity substructure in 
the centre of this galaxy. This galaxy was observed at an 
angle of 46 degrees away from the major axis. 

NGC7012: The radial velocity profile shows a small 
amount of rotation in the centre which might be a KDC. 

NGC7597: The surface brightness profile was published 
by Schombert (1986). The radial kinematic profiles show 
possible evidence for a KDC. 

NGC7647: This is a cD galaxy with a red envelope and 
a large number of small companions (Vitores et al. 1996). 
The surface brightness profile was published by Schombert 
(1986), and HST imaging by Laine et al. (2003). The ra- 
dial velocity profile of this galaxy shows a clear KDC in the 
centre. 

NGC7649: The surface brightness profile was published 
by Schombert (1986), and HST imaging by Laine et al. 
(2003). A KDC is clearly visible in the centre of this galaxy. 

NGC7720: This is a multiple nuclei galaxy and has a 
massive close companion galaxy, also an elliptical. The sur- 
face brightness profile was published by Lauer (1988), and 
HST imaging by Laine et al. (2003). The radial velocity pro- 
file reveals no significant rotation but there might be some 
substructure in this galaxy. 

NGC7768: Surface brightness profiles were published by 
Schombert (1986), Malumuth & Kirshner (1985) and Jordan 
et al. (2004) and HST imaging was published by Laine et 
al. (2003). The radial velocity profile of this galaxy shows 
significant rotation of the order of 114 ± 11 km s _1 along 
the major axis. This agrees with the 101 ± 5 km s _1 and the 
99 km s _1 rotation found previously by Fisher et al. (1995a) 
and Prugniel & Simien (1996), respectively. 

PGC026269: Also known as Hydra A. The radial kine- 
matic profiles of this galaxy reveal rotation of the order of 
51 km s _1 along the major axis, as well as an increasing 
velocity dispersion profile. 

PGC044257: The data indicate that the velocity disper- 
sion might increase with radius, although the large scatter 
and the small number of bins make it difficult to assess the 
slope. 

PGC071807: No significant rotation or substructure can 
be seen in this galaxy. 

PGC072804: The surface brightness profile was pub- 
lished by Malumuth & Kirshner (1985). The radial velocity 
profile of this galaxy shows rotation of the order of 50 km 
s _1 along the major axis. 

UGC02232: This galaxy has an extended envelope and 
a large number of objects in the vicinity. HST imaging was 
published by Laine et al. (2003). No significant rotation or 
clear substructure is detected. 
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UGC05515: The surface brightness profile was pub- 
lished by Schombert (1986). The radial velocity profile shows 
a KDC in the centre of this galaxy. 

UGC10143: This galaxy has a chain of companions. 
The surface brightness profile was published by Schombert 
(1986), and HST imaging by Laine et al. (2003). No signifi- 
cant rotation or substructure can be inferred from the radial 
kinematic profiles of this galaxy. 



6 DISCUSSION 

(i) Five out of 41 BCGs (ESO349-010, ESO444-046, 
ESO552-020, NGCS311, PGC026269) were found to have 
a positive velocity dispersion gradient. 

The radial kinematic studies done so far on early-type 
galaxies are mostly limited to normal ellipticals, for which 
flat or decreasing velocity dispersion profiles are found. The 
majority of the results previously obtained for very small 
samples of BCGs are similar to those of normal ellipticals. 
Fisher et al. (1995a) found one galaxy (IC1101) in their sam- 
ple of 13 BCGs with a positive velocity dispersion gradient. 
Carter et al. (1999) found one (NGC6166) of their sample of 
three BCGs to have a positive velocity dispersion gradient, 
although Fisher et al. (1995a) did not find it for this galaxy. 
Brough et al. (2007) found negative velocity dispersion gra- 
dients in five out of their sample of six brightest cluster and 
group galaxies (the other one had a zero velocity dispersion 
gradient). 

Both IC1101 and NGC6166 form part of the sample of 
BCGs studied here, but in both cases the measured veloc- 
ity dispersion profiles do not reach the radius achieved in 
the much smaller samples in the above mentioned studies. 
Thus, the positive velocity dispersion gradient could not be 
confirmed for IC1101 or NGC6166 (see Figures lA"3l and EH) . 
However, five other BCGs were found to have a positive ve- 
locity dispersion gradient, though admittedly the slope is 
marginally positive in some cases. If these positive velocity 
dispersion gradients are not caused by systematic errors in 
the various data reduction or velocity dispersion measure- 
ments by different authors, then they imply a rising mass- 
to-light ratio. 

(ii) At least 12 BCGs (NGC3842, NGC4889, 
NGC7647, ESO488-027, IC5358, MCG-02-12-039, 
NGC1713, NGC2832, NGC4839, NGC6269, NGC7649 
and UGC05515) show clear velocity substructure in their 
profiles. 

From studies of elliptical galaxies in high density envi- 
ronments (e.g. Koprolin & Zeilinger 2000), the incidence of 
KDCs is observed to be about 33 per cent, rising to 50 per 
cent when projection effects are considered. Hau & Forbes 
(2006) finds KDCs in 40 per cent of their isolated ellipti- 
cal galaxies. For BCGs, we have found at least 12 of the 
41 BCGs show clear velocity substructure, amounting to 29 
per cent of the sample (intermediate and minor axis data 
included) . 

KDCs can be the result of a merger event (Koprolin & 
Zeilinger 2000), but can also occur when the galaxy is tri- 
axial and supports different orbital types in the core and 
main body (Statler 1991). The fact that the incidence of 
KDCs in BCGs compares with that found for normal ellip- 
tical galaxies in high density environments suggests that the 



two classes share the same fraction of galaxies with triaxial 
shapes. 

(iii) NGC6034 and NGC7768 possess significant rota- 
tion (134 an d 114 km s" 1 respectively) along the major 
axis. Several other BCGs show rotation that is > 40 km 
s _1 and more than three times the standard error: ES0346- 
003, GSC555700266 and NGC4839 (major axis spectra); 
ESO488-027, IC5358 and UGC05515 (intermediate axis 
spectra), as do the two elliptical galaxies NGC4946 (major- 
axis spectra) and NGC604-7 (intermediate axis spectra). 

Carter et al. (1999) found small rotation along the ma- 
jor axis at large radii (30 - 40 arcsec) for their sample of 
three BCGs, which is consistent with the nearly complete 
lack of rotation found near the centres of the sample of 13 
BCGs by Fisher et al. (1995a). According to Fisher at al. 
(1995a), the lack of rotation found in samples of BCGs are 
in agreement with the expectation of declining importance 
of rotation with increasing luminosity for elliptical galax- 
ies. The lack of rotation is also compatible with the idea 
that these objects formed through dissipationless mergers 
(Boylan-Kolchin et al. 2006). The remnants left by mergers 
with or without dissipation are expected to differ in their 
kinematical structure. In a merger which involves gas-rich 
galaxies, the gas will form a disk. After the gas has been 
removed from the system at the end of the merger (through 
ejection and converted into stars), the remnant will show 
rotation (Bournaud et al. 2005). Whereas in a merger where 
dissipationless processes dominate, the remnant will show 
little or no rotation (Naab & Burkert 2003; Cox et al. 2006). 
In this study, clear rotation above 100 km s _1 was found for 
NGC6034 and NGC7768, while most BCGs showed little or 
no rotation. This kinematical differentiation (the existence 
of slow and fast rotators) in early-type galaxies is also clearly 
visible in the SAURON data presented by Emsellem et al. 
(2006). 

The amount of flattening that is expected due to rota- 
tion in a galaxy depends on the balance between ordered 
and random motions, and this can be quantified using the 
anisotropy parameter. The rotation of elliptical galaxies is 
conventionally expressed as the aniso tropy p arameter, de- 
fined as (Vmax/cro)* = (V max /ffo) / ye/1 6 (Kormendy 
1982), where the rotational velocity Vm ax is measured as in 
Section O and the central velocity dispersion ctq is taken as 
the measurement for the central velocity dispersion in Table 
[4] A value of (Vmax/co)* ~ 1 would be expected if a galaxy is 
flattened by rotation. The anisotropy parameter (Knax/oo)* 
can be used to separate galaxies that are rotationally sup- 
ported from those that are supported by a anisotropy where 
the value is substantially less than unity. The division oc- 
curs at about (V ma x/cro)* = 0.7 (Bender, Burstein & Faber 
1992). 

Figure [8] shows the anisotropy parameter as a function 
of galaxy B-band luminosity, where the above mentioned di- 
vision is indicated by the horizontal dashed line. Note that 
the errors indicated are only propagated from the errors on 
the velocity measurements taken to be the extreme radial 
velocity points, and the errors on the central velocity disper- 
sion. They do not take into account the general uncertainties 
involved in determining the most extreme velocity measure- 
ments. Therefore, care has to be taken when interpreting 
individual points on the diagram. The only notable BCG 
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Figure 8. The anisotropy-luminosity diagram. The horizontal 
dashed line separates the rotationally supported galaxies from 
the anisotropic galaxies as described in the text. Only the BCGs 
for which major axis spectra were taken (within 10 degrees) are 
plotted. The box plotted in the figure outlines the region contain- 
ing data on giant ellipticals from Bender et al. (1992). 



data point that lies significantly above the dashed line is for 
the galaxy PGC026269, which possesses moderate rotation 
(51 km s _1 ) and surprisingly low central velocity dispersion 
(222 km s -1 ). However, the ellipticity of this galaxy is zero 
and it is not rotationally supported. 

Another factor that complicates the dynamical inter- 
pretation of individual points is that the observed ellipticity 
is a global property of the galaxy, whereas the kinematic 
measurements taken here only reflect the kinematics along 
the axis where the slit was placed, and only close to the 
centre of the galaxy. For example, a disk component may 
dominate the measured kinematics but will have little effect 
on the ellipticity, making the galaxy appear to rotate faster 
than its global ellipticity would suggest (Merrifield 2004). 

For comparison, the sample of isolated ellipticals from 
Hau & Forbes (2006) are also plotted in Figure [8] Of this 
sample, 11 galaxies were observed along the major axis, and 
the central velocity dispersions were derived from the bins 
closest to the galaxy cores. The BCGs show less rotational 
support than the isolated elliptical galaxies, as a class. Spiral 
bulge (typically rotationally supported) and giant elliptical 
data from Bender et al. (1992) are also plotted. Their central 
velocity dispersions were derived over the whole half-light 
radii of the galaxies. Even though large rotation velocities 
were found for a few individual cases, the BCGs are consis- 
tent with the general trend for very massive galaxies. 



7 SUMMARY 

In this paper, we have presented the largest optical spec- 
troscopic sample of BCGs with radial information to date. 
We have derived the rotational velocity and velocity dis- 
persion profiles for 41 BCGs. We find clear rotation curves 
for a number of galaxies for which we have obtained ma- 
jor axis spectra, and two galaxies with rotational velocities 
exceeding 100 km s _1 . However, because of the generally 
large central velocity dispersions, the BCG data are consis- 
tent with the trend for very massive elliptical galaxies on 
the anisotropy-luminosity diagram. At least 29 per cent of 
the galaxies show very clear velocity substructure. 

A future paper will be devoted to the stellar popula- 
tion properties of this galaxy sample. Despite the undeniably 
special nature of BCGs due to their extreme morphological 
properties and locations, the kinematic properties investi- 
gated here (rotation and incidence of velocity substructure) 
seem normal compared with their ordinary giant elliptical 
counterparts. 
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APPENDIX A: KINEMATIC PROFILES: 
COMPARISON WITH LITERATURE 

The kinematic radial profiles of some of the galaxies could be 
compared to previous measurements in the literature. Four 
comparisons are shown in Figures IA1I to IA4I Two of the 
Coma cluster cD galaxies (NGC4839 and NGC4889) were 
compared to data from Fisher et al. (1995a) and Mehlert 
et al. (2000). The third Coma cD (NGC4874) was not com- 
pared as it was observed at a different slit position angle than 
in the other studies. No systematic differences are found 
(within the error bars) for the measurements of this study 
compared to the previous studies. The velocity dispersion 
measurements of NGC4889 in Fisher et al. (1995a) are lower 
than the other studies. This discrepancy was already re- 
ported in Mehlert et al. (2000). NGC6166 was compared to 
Fisher et al. (1995a) and agree within the errors. IC1101 was 
compared to Fisher et al. (1995a), who found the velocity 
dispersion profile to be increasing. As can be seen from Fig- 
ure IA4I the velocity dispersion profile measured here does 
not extend out to the same radius as that from Fisher et 
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Figure Al. NGC4839. Kinematic profile compared to previous 
literature. The data are folded with respect to the centre of the 
galaxy, and the radial velocities are given as relative to the central 
radial velocity. 



al. (1995a) but the central measurements are in good agree- 
ment. 

This paper has been typeset from a TjrjX/ JOTgX file prepared 
by the author. 
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Figure A2. NGC4889. Kinematic profile compared to previous 
literature. 
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Figure A3. NGC6166. Kinematic profile compared to previous 
literature. 



Radial kinematics of BCGs 



R (arcsec) 

2 4 6 8 10 12 

200 i 1 1 1 1 1 1 1 




R (arcsec) 



• This study 

o Fisher et al (1995a) 



Figure A4. IC1101. Kinematic profile compared to previous lit- 
erature. 



